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Abstract 

We propose exclusive difTractive dijet photoproduction as an ideal measure of the ofF- 
diagonal gluon distribution at high scales. We solve the off-diagonal evolution equations for 
the gluon and quark singlet over the full kinematic domain. We discuss the nature of the 
solutions of these equations, which embody both DGLAP and ERBL evolution. We give 
predictions for the transverse momentum distribution of the jets. In particular we quantify 
the enhancement arising from the evolution of the off-diagonal parton distributions. 
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1 Introduction 



Traditionally diffractive processes are described, within perturbative QCD, by two gluon ex- 
change in . An example is the 7p — > qqp process sketched in Fig. |l], where the outgoing qq system 
may emerge either as a vector meson or as two jets. It has been argued that the amplitudes 
for these reactions are proportional to the conventional gluon distribution^ G{x]p) = X]pg{x]p), 
with Xjp given by 

= W^^^ 

where M is the invariant mass of the diffractive qq system and W is the 7p centre-of-mass 
energy. The variable xjp is the fraction of the proton longitudinal momentum transferred to 
the diffractive system by the exchange of the two gluons. A non-zero virtuality of the photon, 
/ 0, allows electroproduction processes to be considered, in addition to photoproduction 
with Q2 = 0. 

We notice that the two exchanged gluons in Fig. |l] carry different momentum fractions x and 
x' , and so the process should actually be described by an off-diagonal0 gluon distribution G{x, x'). 
For vector meson electro- or photo-production it turns out that it is a good approximation to 
use the off-diagonal distribution G(xjp, 0) with x = xjp and x' = 0. The error made in using the 
traditional diagonal distribution G{x]p,xip) = X]pg{x]p) has been quantified in Ref. |5| (related 
work can be found in [Q). There the DGLAP-type evolution equation for the off-diagonal 
distribution G{x,x') was solved and compared with standard DGLAP evolution for xg{x). The 
ratio 

R{x,x') = (2) 
xg{x) 

was calculated as a function of Q^, for different choices of (diagonal) starting distributions. The 
ratio R{x, x') was, as expected, found to be above unity and to increase with and with 
Xjp = X — x' . For example the off-diagonal effects for J /ip photoproduction can be estimated 
from the values obtained for i? at x = xjp and x' = 0. From Ref. |^ we see that the J/ip 
amplitude is enhanced by i?(x,0) ~ 1.1 by off-diagonal evolution. The prediction for the cross 
section is thus increased by a factor of about 1.2. 

Ref. 1^ considered only the off-diagonal effects arising from DGLAP-type evolution and so 
the analysis was restricted to the domain x' >Q. The difference with the conventional diagonal 
distribution was found to be largest when x' = and to decrease rapidly with increasing x' . 
That is ii(x, x') was largest when x' = and rapidly decreased to unity as x' — > x. However the 
DGLAP domain is only a part of a more general evolution of the off-diagonal distribution. In 
the region x' < the gluon G{x,x') looks like the distribution amplitude for the proton to emit 
two gluons and consequently obeys an ERBL-type evolution equation 0, Therefore processes 
which depend on a wider range of x', and in particular sample the x' < region, may be subject 
to much larger off-diagonal effects than the J/V' example that we mentioned above. 

^The scale at which the gluon distribution is evaluated is z{l — z)Q^ +Pt + mg, where and ±pj, are the 
mass and transverse momenta of the emitted quarks, and z is the fraction of photon momentum carried by one 
of the quarks. So perturbative QCD is valid if either (with z ~ 1/2), Px or is large. 

^In the literature alternative nomenclatures are used for "off-diagonal" distributions. In particular non-forward, 
off-forward and non-diagonal are used in Refs. Q, [^] and respectively. These works introduce different, but 
equivalent, definitions of the distributions. Our work uses the formulation of [^, but we prefer to retain the name 
"off-diagonal" since we will be working in the forward direction. 
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Figure 1: The QCD subprocess 7p — > {qq)p', describing the exclusive diffractive photoproduction 
of dijets. 



One such process is the exclusive diffractive production of a pair of (quark) jets with high 
values of transverse momenta ipj-. For this process it is necessary to work in terms of the 
off-diagonal gluon distribution / unintegrated over the transverse momentum of the gluon, 



= (3) 



The computation of the high pT dijet production cross section requires integration over the 
entire region of kx (with an important contribution coming from the region kx ^ Pt)- Hence 
this process explores the detailed properties of the off-diagonal gluon distribution (P) in a broad 
range of x' and kx- In fact we will find that an important contribution comes from the region 
\x'\ ~ xjp which should be compared to x' = which is relevant for diffractive J/V' production. 
As a result we are able to study the much richer structure of the full evolution equations, which 
combine DGLAP-like as well as ERBL-like features. 

The content of the paper is as follows. In Section ^ we present the formula for the diffractive 
dijet cross section driven by the off-diagonal unintegrated gluon distribution. We carefully 
examine the kinematic relations relevant for the off-diagonal analysis. In Section |3| we present 
a simplified analysis which provides valuable insight into the process, as well as forming the 
basis for a comparison with our full off-diagonal treatment. The kinematic relations relevant 
to the perturbative and nonperturbative regions are given is Section Then in Section |5| we 
study in detail the full evolution equations of the off-diagonal distributions, which embody both 
DGLAP and ERBL components. The detailed form of these equations is given in the Appendix. 
The effect of using the off-diagonal distributions to predict the cross section for the exclusive 
diffractive dijet photoproduction is quantified in Section Finally Section contains a brief 
summary of our main results. 

2 General form of the diffractive dijet cross section 

The differential cross sections for the exclusive diffractive production of dijets from transversely 
and longitudinally polarised photons, described by the QCD subprocess j*p — > (qq) p' shown in 
Fig. have the following structure (see also |Q) 



d'^Px dt 



= I I '^!T ^-4^ f{x,x',k^) f{x,x',kl^) <^'^^^{z,kT,kT,PT) , 

(4) 
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where are the transverse momenta of the quark jets, z and (1 — z) are their longitudinal 
momentum fractions, and is the transverse momentum of the gluons. The barred variables 
refer to the amplitude which is the complex conjugate of the amplitude of the process shown in 
Fig. 0. Formula (Q) was obtained assuming that the imaginary part of the amplitude dominates, 
with the neglect of a possible contribution from the real part. To compute the imaginary part 
we have to take into account the four possible ways that the two exchanged gluons can couple 
to the two quarks forming the dijet system. This is reflected both in the structure of the impact 
factors gluon distribution function /. Strictly speaking the notation for the 

arguments of / is symbolic and is only meant to indicate that we are dealing with two gluons 
with different longitudinal momentum fractions x and x' . To be precise / is the following linear 
combination of the off-diagonal distributions @ 



f{x,x',k'j^) — — 



f{x{x'),x',kl) + f{x{x"),x",kl 



(5) 



where the longitudinal momentum fractions {x{x') and x{x") of the first gluon, emitted from 
the proton in Fig. |l|, and x' and x" of the second gluon) depend on whether the second gluon 
couples to the quark with momentum fraction z or (1 — z). For these two configurations we have 
respectively 

_ 4 + 2pT • kT _ 4 + 2pr • kT , . 

~ z{Q^ + " (1 - ^) (Q' + W^) ' 

with 

x(x') = X]p + x' and x{x") = xjp + x" , (7) 
where x^p is specified by (|l]) with a diffractive mass of the dijet system given by 

^2 I ^2 



z{l-z) 



M'' = ' ' l . (8) 



Notice that xip = C,=x — x' (or x — x") plays the role of the asymmetry variable for our process. 
It is also important to note that the momentum fractions of the first gluon are always positive, 
whereas the fractions x' and x" carried by the second gluon may be negative for /cp < 2pT- 
In this case the second gluon is emitted, rather than absorbed as is shown in Fig. |^. For the 
complex conjugate amplitudes we have analogous formulae for x' and x" with k^ replaced by 

kT- 

Finally, the impact factors and are |^] 



and 



where 



$^(z,fcT,'fcT,PT) = [z^ + zf]^i{z,kT,kT,PT)+m^q^2{z,kT,kT,PT) (9) 

^l{z,kTM,PT) = AQ''z''{l-zf ^2{z,kT,kT,PT) , (10) 
i Pt Pt- (Pt + kx) 



^l{z,kT,kT,PT) 



H + Q? H + Q][{PT + kTy + Q 



(11) 



Pt ■ (Pt + kr) _^ {pt + ^t) • {pt + ^t) 



\pI + Ql [{Pt + kT? + Ql [{Pt + kT? + Q'] [{Pt + kT? + Q'] J ' 
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and 



<^2{z,kT,kT,PT) 



+ 



(12) 



[pI + Q] \{Pt + kr? + Q] [{Pt + kr? + Q] [{Pt + kr? + Q 
with 

q'^ = z [1 - z) Q'^ + ml (13) 

Cross section (^) is invariant under the interchange z ^ {1 — z). Thus we may replace the 
upper hmit of the z integration by 1/2 and multiply the final result by 2. Now the upper limit 
corresponds to the minimal value of the pomeron momentum fraction xjpmin (or diffractive mass 
= 4 (pj. + ?Tig)) necessary to produce a dijet system with jets of transverse momenta ^Pt- 
The lower limit Zmin corresponds to the maximal allowed value of xjpmax (usually xjpmax < 0.1). 
Thus from (||) and (|8|) we obtain 




4(p|, + 771,2) 



rCjPmax (Q2 + W^) - Q2 



(14) 



We now study the integrations in (Q) with respect to the azimuthal angles and (j) between 
the jet transverse momentum pj. and the gluon momenta fcy and kx, respectively. The angular 
structure of the impact factors and $2 allows us to rewrite the cross section formula (^) in 
the following remarkably compact forms 

= ^T.4 + - iMz,PT)f + [Mz,PT)f \ (15) 



(Pprpdt 

and 

dcTL 



d^prpdt 



t=o 



f- E^' / . dz{z'{l-zf % \Uz.PT)f\ , (16) 
L Pt J 



where the new "impact factors" (pi and 02 have the following forms 

Mz,Pt)= /#/(x,x-,fc^)l|i^- ^7\ \, (17) 



kj, J '' 2\l + a;2 a + b cos 







Mz,Pt)=I^ fd<t>f{x,x',kl)\-^^- I \. (18) 
kj, J Ll + w^ a + b cos r 







For convenience we have introduced the dimensionless quantities 

uj = Q/ Pt , T = kT/ Pt , a = 1 + + and b = 2t . (19) 



The cross section formulae (15) and ( pig ) are quite general. They do not depend on any particular 
behaviour of the gluon distribution function. The crucial observation to obtain these formulae is 
that the terms which could spoil the new impact factor factorization are odd with respect to the 
angles <j) ^-nd (p. Thus they give a vanishing contribution when the full azimuthal integrations 
are performed. 
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Figure 2: The integrand in the curly brackets in ( |20D as a function oi t = kx/ Pt (sohd curves). 
We also show the contributions from the x' > (dotted curves) and x' < (dashed curves) 
regions separately. The two figures correspond to lo = Q/ pT = ^ (left) and w = 0.1 (right). 



3 Insight from the simpUfied case 



The angular integrations in (|T^) and (|1^) cannot be performed analytically since the arguments x 
and x' of the nonforward gluon distribution / depend on the angle (j). It is illuminating, however, 
to first consider a simplified situation in which / = /{xjp, fc^) depends on xjp = x — x', and not 



on X and x' separately. This situation was implicitly assumed in 
(p and the angular integration can be performed, and we obtain 



Then / is independent of 



and 



Mz^Pt) 



4'2iz,PT) 



IT 



vr 



f{x]P,kT) 



1 n 



—f f{xip,k^) 

Kirp 



62 



(20) 



(21) 



In the limiting case co = Q/ pT ^ 0, which corresponds to diffractive dijet photoproduction with 
high values of pT, the only appreciable contribution to the cross section (|T^) comes from the 0i 



term. In this limit (2 



(1 - T^) and x/(a2 - fe^) 
dk^ 



1 1 — r I and so 



4>i{z,Pt)^-^ J -j^ f{x]P,kT) &{kT-pT)- 



becomes simply 
(22) 



Due to the 1/ k^ factor, the dominant contribution to the integral comes from the region kx ^ Pt- 
In other words the diffractive production of dijets with transverse momentum p^ measures the 
gluons forming the QCD pomeron at kx ^ Pt @ • For kx < Pt the exchanged gluon system is 
unable to resolve the q and q separately which results in zero net coupling and is reflected by 



the theta function in (22). 

In order to better understand why the contribution for gluon momenta /ct < Pt vanishes, 
we divide the integration in (^) into the part corresponding to x' > and the part with x' < 0. 
We find from (^ that the condition x' > is fulfilled if the azimuthal angle (j) between the jet 
prp and the gluon fc^ is given by 



< (b < a 



arccos 



kT 
'2^ 



(23) 
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Figure 3: The cross section da/dp\. as a function of the ]ei pT- The continuous curve shows the 
cross section with both the regions x' > and x' < included. The dashed curve corresponds to 
the contribution from the x' > region only, computed with the infrared cut-off /cto = 1 GeV. 



whereas the region a < 4> < tt corresponds to x' < 0. Notice that if kj< > 2px then only the 



region x' > contributes. For kT < 2pT we split the angular integration in (17) into the positive 
and negative x' parts and carry out the integrations separately 



krp 



+ 



IT 

2 



(24) 



to 



62 



a/3 
vr 



+ 



1 



■ UJ 



l+a;2 



1-^ 
vr 



62 



TT 



where the function /? is defined by 



1 



fh-\- a cos a 

p = — arccos ; 

a \a + b cos a 



(25) 



Notice that the sum of the two contributions in the squared brackets, coming from the x' > and 
x' < regions, simplify in such a way that formula ( pO| ) is recovered. In the limiting case when 
a; = 0, the positive and negative x' parts are equal to a (1 — /?)/7r and —a (1 — /?)/7r respectively, 
and so their sum gives a vanishing contribution to (pi for kx < Pt- As a consequence there 
is no necessity for an infrared cut-off kxQ on the integration in (^). For w 7^ 0, however, the 
cancellation occurs only between the contributions from the two x' regions that are linear in kT- 
The first nonzero term is proportional to u'^ kj^ and leads to a mild logarithmic dependence of 
(pi on kxo- 

The results of the above analysis are well illustrated by Fig. |2| where the solid lines show the 
integrand in the curly bracket of ( pO|) and the dotted and dashed lines show, respectively, the 
x' > and x' < components in (p^). Both the cancellation for r < 1 and the absence of the 
x' < component for kx > 2pT are evident. We take = and px = 10 GeV. The two plots 
correspond to a; = Q/ px = and to a; = 0.1 {rrig = 1 GeV). 

It is informative to show the cross section da/dpx obtained from (Q) for the diffractive 
photoproduction of dijets in the simplified case when the gluon distribution / = 1 and as is 
fixed. The prediction is given by the continuous curve in Fig. g At large px we anticipate from 



I) that 



1/ Px and so da/dp^ ~ 1/px- Such behaviour is evident in Fig. y. Crucial to 
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this behaviour is the cancellation between the x' > and x' < contributions that are linear 
in kT- This is well illustrated by the dashed curve in Fig. ^, which is the cross section which 
would have resulted if only the x' > contribution were included. Then ~ 1/ (pt^to) where 
fcro is the infrared cut-off imposed on the kr integration. As a result the cross section behaves 
as 1/ (Pt^^to) ^^'^ becomes strongly dependent on the value chosen for the cut-off momentum 
kTo- For interest we have extended the predictions in Fig. ^ down to small values of pT- The 
large pT asymptotics which we have described above become apparent for pt ^ 4 GeV. 

So far our discussion has been based on the assumption that the gluon distribution in formu- 
lae (ITtI ) and ( |l8|) does not depend on the azimuthal angle. In this case the both x' regions are 
weighted by the same value of the gluon distribution /(xjp) when the angular integration in ([TtI ) 
and (^) is performed. As a result we have the delicate cancellation in the region of kx < Pt, 
which we emphasized above. This could be significantly changed when the true off-diagonal 
gluon distribution f{x,x') is used, and the two x' regions are weighted differently during the 
angular integration. The study of the influence of this effect on the diffractive dijet cross section 
is a major objective of our analysis. 



4 Pert ur bat ive and nonperturbative contributions 

To evaluate the cross sections given in (jTsI ) and (1T6| ) we have to compute the impact factors (pi 
and (j)2 which depend on the off-diagonal gluon distribution. We divide the kx integration in 
( p^ ) and ( p^ ) into a nonperturbative and perturbative region according to whether kx is smaller 
or greater than kxo ~ 1 GeV. The angular integration can be performed analytically in the 
nonperturbative region. For this purpose we analyse the angular dependence of the arguments 
X, x' and x" of the gluon distribution function (P) by dividing (^) by xjp given by (||) and (|8|) 

x' ,^ , /r2 + 2Tcos0\ x" f t'^ + 2t cos 

= (1 - ^ 7—^ ' = ^ 7—^ • (26) 

X]p V / X]p \ 1 + J 

In the nonperturbative region r = kx/pT 1, and thus we have 

|x'|, \x"\ <C xjp and x « xp . (27) 

These relations allows us to expand the gluon distribution (^ around the point x = xjp, i.e. we 
study / as function of x when the value of xjp is fixed. Thus we obtain 

f{x,x',kl) ^ f{xjp,0,4) + ^^If (^^'O'^t) • (28) 



We also expand the expressions in the curly brackets in (17) and ([Iq ) in powers of r. Then 
the angular integration can be performed analytically in the nonperturbative region. The term 
linear in r vanishes and as a result the first nonzero term is proportional to (that is kj^). This 
allows us to express the resulting kx integration in terms of the off-diagonal integrated gluon 
distribution at scale /c^q using 

/ ^/(xiP,0,4) = G(xjp,0,4o), (29) 
JO k^ 

tflfj^-"-'^^ = f(^r.0.4„)^ (30) 
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Note that the ofF-diagonal gluon distribution G{x,0, k"^) is not to be identified with the con- 
ventional symmetric gluon distribution G{x,x,k'^) for which the longitudinal fractions are equal 
and hence the asymmetry variable xjp = 0. 

Thus we finally obtain the following decompositions of the impact factors (pi and (j)2 into the 
following nonperturbative and perturbative contributions 



+ 



°° dk^ 
fc2g k^ 



d(j) f{x,x',kl) 



a + 6 cos I 



(31) 



and 



vr 



Pt 



UJ 



{I +00 



2^3 



G(xjp,0,4o) + 



1 dG , ,2 , 



~ 2 r 1_ (32) 

For difFractive dijet production with large pT by far the dominant contribution comes from the 
perturbative region. 

Relations (|2^) are not generally valid in the perturbative region. For example for gluon 
transverse momenta kx ^ Pt, which dominate in the simplified case of Section we find that 
at least one of the fractions \x'\, \x"\ ~ xjp. For z = 1/2 and a; <C 1 this is equivalent to 

1/2 X]p < x < 5/2 xjp . (33) 

Thus in the perturbative region we study a much broader range of x of the off-diagonal gluon 
distributions (|5|) with the second gluon momentum fractions x' or x" being both positive and 
negative. Also the ofF-diagonal distributions are required over quite a range of k^ and are 
computed using evolution equations. We discuss the off-diagonal evolution in detail in the next 
section. 



5 Off-diagonal parton distributions and their evolution 

Before discussing the evolution equation for the off-diagonal parton distribution it is convenient 
to follow Radyushkin |lO| and use the notation^ 

G,^{x,kT) = G{x,x',k^) , (34) 

for the off-diagonal gluon distribution, where the dependence on x' is implicit through the 
relation x' = x — xjp. Both the variables x and xjp lie in the interval to 1. Therefore we view 
the off-diagonal distributions as a family of functions (of x and the scale kx) labelled by the 
asymmetry variable xjp. Radyushkin uses C for xjp, but for the moment we will keep the latter 

* Alternative formulations of the off-diagonal evolution equations can be found in M and |L1|. 
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notation to indicate that xjp defined in (||) is the asymmetry variable in the process that we 
consider. 

In the region of x > xjp (that is x' > 0) the ofi'-diagonal gluon distribution Gxp{x) describes 
a gluon emitted by the proton with momentum fraction x together with a gluon absorbed with 
fraction x' = x—xjp. Thus Gxp{x) is a generalization of the ordinary diagonal gluon distribution, 
that is in the diagonal limit Xjp = x — x'^Owe have 

Gxp=o{x) = G{x,x) = xg{x) . (35) 

We stress once again that we must distinguish between G{x,x) and Gx{x) = G{x,0). In the 
first case the asymmetry variable xjp = while in the later xjp = x. 

In the region x < xp the second gluon has a negative fraction x' of the proton's momentum. 
In this case it is more transparent to use —k' for the momentum of the second gluon, that is we 
make the replacement 

k' = x'p + kT k' = i-k') = \x'\p-kT . (36) 

Thus we have the two gluons emitted form the proton with the momenta k = xp + kx and k' . 
The gluon distribution Gxp{x) for x < xjp may be regarded as the probability amplitude for 
emitting a two gluon colour singlet system from the proton with the individual gluons carrying 
fractions x/ xp and \x'\/ xp of the total momentum of the system r = xpp = k + k' . 

To evaluate the perturbative contribution we use the full form of the off-diagonal evolution 
equations for the gluon distribution Gxp{x,fi) {fi = kx in our case) and the off-diagonal singlet 
quark distribution 

Nj 
i=l 

where g^.^ and qt^^ are off-diagonal quark and antiquark distributions introduced in analogy to 
the diagonal case. The evolution equations have the following general form 



^xpix,fi) = J dz P^^{x,z;fi) T.xpiz,n) + j dz P^^{x,z;n) Gxpiz,^) 



d_ 



(38) 



Gxp{x,n) = J dz P^^{x,z;fi) T.xp{z,n)+ J dz P^^ {x , z; fi) Gxp{z,fi), 



where the structure of the kernels P depends on the relations between the x, z and xp variables. 
The detailed form of the evolution equations extracted from the review is given in the 
Appendix. 

For X > Xp {x' > 0) the integration in ( p8| ) covers the range from x to 1 and the kernels have 
a form which coincides with the standard Altarelli-Parisi kernels for xp = 0. That is (|3^ ) reduce 
to the DGLAP evolution equations [^] in the limit xp = 0. On the other hand for x < xp 
[x' < 0) the form of the kernels depends on whether the integration variable z is less or greater 
than X. In this case the integration covers the whole range (0, 1) and in the the limit xp = 1 
Eqs. (|3^) reduce to the ERBL evolution equations 0, for the parton distribution amplitudes 
in a meson. 
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Figure 4: The off-diagonal singlet and gluon distributions and the logarithmic derivative of the 
gluon distribution (dashed curves) at = 10^ GeV^ for C, = xp = 10~^ (middle vertical line). 
The continuous curves show the initial distributions at /Xq = 1 GeV^ (lower curves) and effect 
of their evolution using the DGLAP equations (upper curves). 

We solve the evolution equations numerically using an expansion in terms of Chebyshev 
polynomials and cross-check the solutions against, on the one hand, the known asymptotic 
solutions of the ERBL equations and, on the other hand, the solutions of the DGLAP evolution 
equations. Moreover particular attention is paid to the behaviour of the resulting off-diagonal 
parton distributions at point x = xjp, where the two (DGLAP-like and ERBL-like) forms of 
the off-diagonal evolution equations coincide. In all cases we find excellent agreement with the 
expectations resulting from analytical insight. 

In the two upper plots of Fig. ^ we show the off-diagonal singlet and gluon distributions 
evolved up to fi^ = 100 GeV^ (dashed lines) starting from initial conditions specified at fiQ = 
1 GeV^ (lower solid curves). We postulate the following form of the initial distributions 

^xipix) = (1 - xjp)" { (1 - xip) SMRs(a;) + xp T.bl{x) } , 

(39) 
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Gxpix) = {1 - xip)"' { {1 - xip) Gmrs{x) + xp Gbl{x) } , 

where Smrs aiid Cmrs the MRS distributions obtained from global fits based on the 
DGLAP evolution equations to data for deep inelastic and related hard scattering processes, 
and Ebl ~ a; (1 — x) and Gbl ~ (1 — x)^ are asymptotic solutions of the ERBL equations (for 
the nonsinglet and pure gluon cases). The postulated form in the curly brackets encompasses 
the basic feature of the off-diagonal parton distributions (and evolution equations) that in the 
limits xjp — > or 1 the DGLAP or ERBL components, respectively, are obtained. However in 
case of a nucleon it is unlikely that the two partons can be emitted, sharing the whole nucleon 
longitudinal momentum (which corresponds to the xjp = 1 limit). This observation is accounted 
for by the additional powers of (1 — xjp) with n,m > 0. The asymmetry variable was chosen 
to be xjp = 10~^, and is indicated by the middle vertical dotted line (the outermost vertical 
lines correspond to the values xjp/2 and 2xjp). Clearly the initial distributions are dominated 
by the DGLAP (MRS) component for such a small value of xjp. We may compare the dashed 
curves (the evolved off-diagonal parton distributions) with the upper continuous curves which are 
obtained from the same initial conditions using the conventional DGLAP evolution equations^. 
The bottom plot compares the logarithmic derivative of Gxp{x), which is used in our analysis, 
see (^), at the same scale //^ = 100 GeV^. Again the dashed curve corresponds to the off-diagonal 
case while the continuous curve is obtained using the DGLAP evolution equations. 

The solutions reflect the mixed (DGLAP and ERBL) nature of the off-diagonal evolution. 
This is particularly visible for E(^=a;^(x) which becomes negative for x < xjp. This feature is 
found for the solution of the ERBL evolution equations for the distribution amplitude which is 
obtained from S^(x) in the limit C = xjp = 1. Thus the off-diagonal parton distributions do not 
have a probabilistic interpretation, unlike the conventional diagonal distributions. The region 
|x — xjp| ~ xjp, indicated by the three vertical lines, is particularly interesting from the point of 
view of the difference between the off-diagonal and diagonal parton distributions. This region is 
relevant for gluons with momenta kx ~ Pt (see (|33|)) which are supposed to give an important 
contribution to the dijet production cross section. We will explore this effect in the next section. 



6 Effect of off-diagonal distributions on the dijet cross section 

For diffractive dijet photoproduction only the cross section ( [T5| ) for the transversely polarised 
photons contributes. In addition the (p2 part is strongly suppressed for high values of pr- Thus 
we need only analyse the impact factor given by (^T]). In terms of r = kx/pT the perturbative 
part of (pi has the following form for u = 



^^r'\^,PT) 



27r 
Pt 



dr 



TO 



JO 



— fxp{x,TpT) -< 1 
TT Z 



1 + + 2t cos ( 



(40) 



where tq = kxo/PT- The function fxp is written in the notation introduced by Eq. (^) and is 
related to the off-diagonal gluon distribution Gxp through relations (^) and (|5|) . In Fig. |5| we 
show the integrand in the squared brackets in ( |40| ) as a function of r. We choose for illustration 
Pt = 5 and 10 GeV for dijet transverse momentum. For each choice we assume the minimum 
value of xjp = 4pj,/ W'^, which corresponds to z = 1/2. 

^The off-diagonal gluon distribution (Fig. ^) in the a;' > region, to the right of the central vertical line, 
together with the conventional DGLAP curve reproduce the main features of the earlier study of Ref. 
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Figure 5: The integrand in the square brackets of (40) for two values of dijet momentum pT 
and the corresponding minimal values of xjp. The continuous curves show the integrand in the 
simplified case when f^p = fxpixjp), while the dashed curves show the full off-diagonal case 
fxp = fxp{x)- The dotted curves adjacent to the dashed ones show the approximate integrand 
of (41) when only the term linear in r is retained. 



The results depend subtlely on the properties of the gluon distribution f^p- Let us start 
with the simplified case discussed in Section]^. That is we assume that fxp = fxpixip^TpT) in 
(^), hence the gluon distribution does not depend on the azimuthal angle. DGLAP evolution 
is appropriate in this case and the resulting integrand in ( pO| ) is shown by the continuous curves 
in Fig. ^. The angular integration leads to the step-like form of the integrand which is then 
modified by the 1/r^ factor to give a function peaked at r ~ 1. This result was anticipated in 
Section ^. 

Now let us explore what happens to the integrand when we use the full off-diagonal gluon 
distribution fxpix^rpx), evolved using Eqs.(|3^ from the input given in (^9|). The results 
are shown by the dashed curves and are strikingly different from the continuous curves of the 
simplified case. First we see that for r > 1 the integrand is larger. This can be easily understood 
by inspecting the off-diagonal gluon distribution shown in Fig. If r > 1 (that is kx > Pt) 
then during the angular integration the argument x of the gluon fxp{x) spans approximately 
the region indicated by the outer vertical lines with the gluon following the dashed curves, see 
also (^). In contrast in the simplified case, the function fxp{x]p) is constant as x spans the 
allowed range, with its value given by where the solid curve crosses the central vertical line in 
Fig. The difference in these forms of fxp explains why the dashed curve in Fig. ^ is larger 
than the continuous curve for r ^ 1. 

The biggest change is in the region r < 1. Recall that in the simplified case the regions 
of x' > and x' < (to the right and to the left, respectively, of the central vertical line in 
Fig. ^) are weighted by the same value fxp{x]p) when the angular integration is performed. As 
a result the terms linear in r coming from the two x' regions cancel (see Fig. ^ for illustration 
of this mechanism). However the off-diagonal distributions are asymmetric with respect to the 
line X = xp which divides the two x' regions (see dashed curve in Fig. ^). This leads the lack 
of the above mentioned cancellation and in consequence to an important contribution for r < 1, 
shown by the dashed lines in Fig. |5[ We quantify this observation by expanding the function in 
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Figure 6: The difFractive dijet photoproduction cross section (in pb/GeV^) integrated over the 
interval xjpmin < xjp < 0.05 as a function of the jet transverse momentum (in GeV). The upper 
triangular points are obtained in the off-diagonal analysis and the lower ones correspond to the 
simplified case. The open points are obtained for fxpix) evolved with the DGLAP equations. 



the curly brackets in (M) in powers of r 



^(pcrt) 



27r 



(r < 1) ~ ^ I dr 



Pt 



fxp{x,TpT) {tCOS< 



cos 20 + ... } 



(41) 



TO 



Now we retain the term linear in r and perform the angular integration. The resulting integrands 
are the dotted lines in Fig. ^, extended into the r > 1 region to make them visible. Thus in 
the off-diagonal case the linear in r term from the above expansion does not vanish after the 
angular integration as it did in the simplified case. Moreover the whole contribution from the 
T < 1 region is essentially determined by this term. 



At this point we may worry that the contribution coming from the linear in r term in (41) 
drastically increases the sensitivity of cj)^^^^^ , and in consequence the dijet cross section, to the 
choice of the infrared parameter kxa- Fortunately this does not happen. The reason is as follows. 



For To < r <C 1 we may additionally expand the function fxp{x) in (41) around x = xp using 



formula (28). Then we obtain 



27r r dr 
Pt 



f dr f d(j) ( dfxjp{x,TpT) 



TO 



2tt 



dr 



1 dfxjp{x,TpT) 



Px J TO 2r 



din: 



X=X]p 



X=Xjp 



T COS ( 



(42) 



where the angular integration has been performed to give the final result. Notice that in the 
simplified case of Section ^ only the first term in the expansion of fxp exists and the vanishing 
contribution is obtained after the angular integration. The full off-diagonal treatment introduces 
an additional linear in r term which leads to the first nonvanishing contribution being quadratic 
in r. As a result (p^^^^^ still depends at most logarithmically on the parameter kxa- We have 



13 



checked that the integrand in (^) reproduces to a good approximation the exact integrand in 
(^ ) up to T ~ 0.3 for both the values of pT chosen for Fig. ^. 

The above analysis also shows that the size of the negative contribution to (f)^^'^'^^^ is mostly 
determined by the value of logarithmic derivative in x of fxjp{x), taken at the point x = xjp. In 
Fig. ^ this derivative corresponds to the slope of the dashed curve at the intersection point with 
the middle vertical line. On inspecting Fig. |^ notice that the gluon distribution evolved with 
the DGLAP equations (continuous curve) is also asymmetric with respect to the line x = xjp. 
Thus, when such gluon is used in the dijet analysis, a similar negative contribution will be 
present. The size of this contribution, however, is smaller than in the fully off-diagonal case 
since the slope of the continuous curve at x = xjp is always smaller than that of the dashed 
(off-diagonal) one. 

In summary, the off-diagonal parton distributions enhance the contribution to (pi in com- 
parison to the simplified case for > Pt and lead to a significant negative contribution in the 
region kx < Pt- Both contributions reflect the particular form of the gluon distribution (|3|) that 
is dictated by the off-diagonal evolution equations (|3^). We have checked that for pt > 5 GeV 
the net effect after the r (fcy) integration in (^) is always negative in the whole range of the 
z variable (zmin < z < 1/2) in the fully off-diagonal case. Moreover, after squaring and inte- 
grating over z in (p^, the impact factor (pi computed using the true off-diagonal distributions 
leads to a significantly larger cross section for diffractive dijet photoproduction ([l5| ) than in the 
simplified case. In Fig. ^ we show the cross section da/dpqp, integrated over xjp between the 
values xjpmin = ^p\/W'^ and xjpmax = 0.05. In addition we approximate the effect of the inte- 
gration over the momentum transfer t by dividing ( [l5| ) by the diffractive slope A^i = Q GeV~^. 
The upper triangular points correspond to the true off-diagonal result while the lower ones are 
obtained in the simplified analysis. The open points show the result when fxjp{x, kx) is evolved 
using the DGLAP equations. 

We emphasize the special nature of the exclusive diffractive photoproduction of dijets as 
a probe of the off-diagonal gluon distribution. In our study we have assumed that the cross 
section is dominated by the QCD subprocess 7p iQQ)p'- A possible contamination may arise 
from the subprocess 'jp iQQg)p', in which the qq pair and the gluon form the two jet system. 
This background process has been estimated in |14| for the inclusive diffractive cross section and 



found to be important in the region of low pT and large diffractive mass M, but its contribution 
should be small in our case since we are studying the large px domain. 



7 Conclusions 

The diffractive photoproduction of dijets with high pT is an ideal probe of the properties of the 
off-diagonal (unintegrated) gluon distribution f{x,x',k^) over a wide kinematic range. Indeed 
the calculation of the differential cross section da/dpj^ involves the integration over entire range 
of kx with important contributions coming from the region kx ^ Pt, and so requires the solutions 
of the off-diagonal evolution equations to determine /. Our detailed studies of these equations 
revealed some novel features. We found that the region kx < pt gives an important contribution 
to the cross section, originating from the specific properties of the evolved gluon in the x' < 
(ERBL) and x' > (DGLAP) domains. This is in contrast to the simplified treatment, discussed 
in Section ^, in which the entire contribution to the cross section comes from the region kT ^ Pt 
leading to a far smaller cross section. In the simplified analysis the gluon distribution was 
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assumed to be only a function of xjp (and kj^) - an equivalent approximation was made in the 
analysis of Ref . fH . 

To summarize we found that if the true off-diagonal gluon distribution is used then the impact 
factor (Eq. (|3l|)) receives an important negative contribution from the /ct < Pt region as 
well as an enhanced positive contribution from the kx > px domain. The negative contribution 
increasingly dominates with growing pT- After integration over kT the net effect is an overall 
enhancement of the dijet cross section in comparison with either the simplified case or with the 
case in which the conventional DGLAP evolution is used. For example, the enhancement of 
da/dp'^ at pT = 10 GeV due to off-diagonal effects is a factor of 3 as compared to using DGLAP 
evolution, and a factor of 6 compared to the simplified treatment of Section ^ and Ref. Q. The 
corresponding factors for jets of px = 5 GeV are 1.5 and 2.2 respectively. These enhancements 
shown in Fig. ^. can, in principle, be tested experimentally offering an ideal testing ground of 
off-diagonal effects. 
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Appendix 



Here we present for reference the full form of the evolution equations for the off-diagonal singlet 
and gluon distributions and following the prescriptions for the kernels given 

by Radyushkin in [^]. We use his definitions and notations of off-diagonal (asymmetric in his 
language) parton distributions. In this notation is the asymmetry (skewedness) parameter 
which in the main body of the paper is equal to the variable xjp 



C, = X — x' = xp , 



(43) 



with < C < 1. The distributions are defined in such way that in the limit C ~^ (the DGLAP 
limit) we obtain the ordinary parton distributions 



S^(x,/x) 



G^(x,/i) 



{q'{x,fi) +q\x,n)} 



i=l 



X gix,n) . 



(44) 



(45) 



Notice that the ordinary quark distributions are not multiplied by x as in the gluon case. In the 
opposite limit (" ^ 1 (the ERBL limit) the off-diagonal distributions become the distribution 
amplitudes for finding two partons in a meson sharing fractions x and 1 — x of its momentum. 

^^{x,^^) ^ ^ {vI/Xx,;u)+*\x,^)} (46) 



1=1 



Gi;{x,n) 



(47) 



Using the notation x' = x — C and z' = z — we have the following equations for x > ^ (that is 
for x' > 0) 



d 



dz 



X — z 



f + 7) S,(x,/.)-(l + ^)S,(.,/.) 



3 ^ (1-X)2 

- +\n- — 

2 1 - C 



+ ^^Nf j — 
vr J zz' 

X 



z J \ z' J zz' 



G(;{z,n) , 



(48) 



d 



fj.-^Gr{x,fj.) = ^^^^^ Cf I dz 

Of! IT J 



T,(;{z,n) 



+ 



TT 



+ 



[{x/z) + {x'/z')] G^x,/,) - [{x/zy + {x'/z'y] Gc{z,^,) 



+ Gi;{x,fi) 



X — z 

11 - 2/{3Nf) 

2Nr 



+ In 



il-x) 



1-C 



(49) 
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where Cp = 4/3 and = 3, and A'^^ is the number of active flavours. In the hmit C = the 
above equations become the DGLAP evolution equations. 

The equations for x < ^ (that is for x' < 0) are more complicated since they involve 
integration with different kernels in the intervals (0,x) and (x, 1). We have 



vr 



Cf { dz 



C 



X — z 



dz 



C 



+ 



z — X 



3 , x(l — x) 

- +ln— 

2 C 



Nf 



IT 



dz / x' 



dz f x 



X 2x — C 
4- + - ^ 

C C- z 



4(1-^1 + 



C-2x 



Gdz,^^) 



(50) 



M^G'c(^'/^) 







1-:^) ^dz,fi) + I dz(^2-^) ^dz,fi) 



+ 



as{^l) 



Nr 



2 / x\ f X X 



Gdz,ij) 



+ 



{x'/z') Gdz,fi)-Gdx,fi) 



X — z 



1 

+ J dz 

X 



2x ^3_2^ + C^^ 



zjiev C 

+ Gdx,ij) 



z — X 



ll-2/(3iV^) ^^^x(l-x) 



C 



(51) 



For C = 1 the above equations reduce to the ERBL evolution equations for the distribution 



amplitudes. It is also instructive to check that both set of equations, that is (48,49) and (5C,51), 
lead to the same limiting set of equations when x — > C from both sides. 

The four equations (48-51) form a coupled set of equations which, in general, need to be 
solved simultaneously. However for x > C it is enough to solve the first two since the right hand 
side of these equations involves parton distribution for values z > x (as is true for the DGLAP 
equations in the limit C = 0). This is not the case if x < C,. The solution depends on values of 
parton distributions in the full interval (0, 1) and we have to solve all four equations together. 
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